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Abstract 

We analyze the production of the neutral Higgs particles of the Minimal Supersymmet- 
ric extension of the Standard Model at the Fermilab Tevatron collider. We consider the 
two main production and detection channels: gluon-gluon and bottom quark fusion lead- 
ing to Higgs bosons which subsequently decay into tau leptons, gg, bb ^Kiggs ^ t^t~ . 
We update the production cross sections and the decay branching ratios and obtain pro- 
duction rates that are significantly smaller at high masses than the ones used by the 
CDF and DO experiments in their search. We then evaluate the various theoretical un- 
certainties that affect these rates, uncertainties that have not been been considered in 
the CDF/DO analyses and which turn out to be rather large. Including these two effects 
will significantly loosen the constraints obtained on the supersymmetric Higgs sector at 
the Tevatron. 



The search for the Higgs bosons, the remnants of the spontaneous breaking of the elec- 
troweak symmetry that is at the origin of the elementary particle masses, is the main goal of 
present high-energy colliders. While a single Higgs boson is predicted in the Standard Model 
(SM), the minimal realization of the symmetry breaking with only one Higgs doublet field [T], 
the Higgs sector is extended in supersymmetric theories [2], that are widely considered to be the 
most attractive extensions of the SM as they stabilize the hierarchy between the electroweak 
and Planck scales induced by the large radiative corrections to the Higgs boson mass. In the 
minimal extension, the Minimal Supersymmetric Standard Model (MSSM) [2], two Higgs dou- 
blet fields are required, leading to the existence of five Higgs particles: two CP-even h and iJ, 
a CP-odd A and two charged particles [3111]. 

In the MSSM, two parameters are needed to describe the Higgs sector at tree-level: the 
mass Ma of the pseudoscalar boson and the ratio of vacuum expectation values of the two 
Higgs fields, tan/3, that is expected to lie in the range 1 < tan/3 < 50. At high tan/3 values, 
tan 13 > 10, one of the neutral CP-even states has almost exactly the properties of the SM Higgs 
particle: its couplings to fermions and gauge bosons are the same, but its mass is restricted 
to values MJ^^^ ^ 110-135 GeV depending on some SUSY parameters that enter the radiative 
corrections [1] . The other CP-even and the CP-odd states, that we will denote collectively by 
$ = ^4, H{h), are then almost degenerate in mass and have the same properties: no couplings to 
gauge bosons, while the couphngs to isospin down-type (up-type) quarks and charged leptons 
are (inversely) proportional to tan/3. 

Thus, for tan/3> 10, the $ boson couplings to bottom quarks and r-leptons are strongly 
enhanced while those to top quarks are suppressed. As a result, the phenomenology of these 
states becomes rather simple. To a very good approximation, the $ bosons decay almost 
exclusively into bb and t~^t~ pairs with branching ratios of, respectively, ~ 90% and ~ 10%, 
while the other decay channels are suppressed to a negligible level [5]. The main production 
mechanisms for these particles are those processes which involve the couplings to bottom 
quarks. At hadron colliders, these are the gluon-gluon fusion mechanism, gg ^ ^, which 
dominantly proceeds through 6-quark triangular loops [6l[7] and bottom-quark fusion, bb ^ ^ 
[814T0]. in which the bottom quarks are directly taken from the protons in a five active flavor 
scheme [TT]. The latter process is similar to the channel pp bb^ when no 6-quarks are 
detected in the final state |llj . 

With its successful operation in the last years, the Fermilab Tevatron collider has now 
collected a substantial amount of data which allows the CDF and DO experiments to be sensitive 
to the MSSM Higgs sector. Stringent constraints beyond the well established LEP bounds 
M^,M/j > Mz and tan/3 > 3 [12], have been set on the MSSM parameter space [iifA5tan/3] 
using the process gg, 66 ^ $ — t- r+r~. Moderate A masses. Ma ~ 100-200 GeV, together with 
high tan/3 values, tan/3>30, have been excluded at the 95% confidence level (CL) [T3|[T^. 

Nevertheless, a very important issue has been overlooked in these experimental analyses: 
the theoretical uncertainties that affect the production and decay rates, which can be important 
despite of the fact that some higher order perturbative corrections to these processes are known. 
These are mainly due to the unknown higher order corrections in perturbation theory, the still 
not satisfactory parametrization of the parton distribution functions (PDFs), as well as the 
parametric uncertainties stemming from the not very precisely measured values of the strong 
coupling constant as and the bottom quark mass iVffe. In a recent analysis [15], it has been 
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shown that these uncertainties can be rather large at the early stage of the CERN large 
Hadron Collider (IHC), in much the same way as in the case of SM Higgs production both at 
the IHC [15] and at the Tevatron [16]. 

In this Letter, we first update the decay branching ratios and production cross sections 
of the $ bosons and find that the latter are significantly lower at high masses than the ones 
assumed in the Tevatron analyses. We then evaluate the theoretical uncertainties that affect 
these rates and find them to be very large, possibly lowering the cross sections times branching 
ratios by a factor 2. When included in the DO and CDF combined analysis of the MSSM 
Higgs bosons [H], the correct normalization and the uncertainties will drastically reduce the 
MSSM [M^,tan/3] parameter space that has been excluded. 

For the evaluation of the cross sections in the gg ^ ^ and hh ^ ^ production processes at 
the Tevatron, we will concentrate on the pseudoscalar A case and follow very closely the recent 
analysis performed in Ref. [15] for the IHC. We calculate (T{gg — )■ A), known up next-to-leading 
order (NLO) only [7], using the program HIGLU [T7] with central values for the renormalization 
and factorization scales, fiR = fip = fiQ = \Ma, only the dominant loop contribution of the 
bottom quark is taken into account. For the bb ^ A process, known up to next-to-next-to- 
leading order (NNLO) [10], we use the program bbh@nnlc0 with a central scale Hr= jJip = 
fiQ = \Ma- In both cases, we work in the MS scheme for the renormalization of the 6-quark 
mass; however, while rnhifrih) is used in the gg process, frih^fiR) is adopted in the bb channel. 
The resulting partonic cross sections are then folded with the latest MSTW sets of PDFs [IB] , 
consistently at the respective orders, NLO or NNLO, in perturbation theory. 

In both the gg ^ A and bb ^ A processes, we assume the bbA coupling to be SM-like, 
^Abb = TTLh/v. To obtain the true cross sections, one has to rescale the numbers which will 
be given by a factor tan^ /3. In addition, to obtain the cross section for both the A and H{h) 
bosons, an additional factor of two has to be included. As a consequence of chiral symmetry for 
M$^m5 and since the H{h) masses and couplings are very close to those of A, this represents 
an excellent approximatioiu. 

The results for the cross sections a{gg^A) and a{bb^A) are shown in the main frames of 
Fig. 1 for the Higgs mass range that is relevant at the Tevatron, Ma = 90-200 GeV. We have 
compared our values with those given by the program that has been used by the CDF and DO 
collaborations for their cross section normalization, FeynHiggs [22]. This program, initially 
supposed to only provide precise values for the MSSM Higgs masses and couplings, gives also 
grids for production cross sections which should be used with care. For the bb^A channel, we 
obtain cross sections that are ~ 30% smaller. The reason is that FeynHiggs simply provides 
the values given in the original paper [10] which uses the outdated MRST2002 set of PDFs 
which are only partly at NNLO. In the case of gg — ?■ A, the agreement is better as we obtain 

^We thank R. Harlander for providing us with his code. 

^Note that there are additional SUSY contributions in gg — >■ H/ h which do not appear in gg-^ A but they 
are very smah for a large SUSY breaking scale, M$ ^ Ms [E]. Furthermore, there are one-loop vertex 
corrections to the $&6 coupling due to SUSY particles which can be significant as they grow with tan/3 [20] . 
They are implemented in the major codes which calculate the MSSM Higgs spectra and can be readily included. 
However, in the case of pp— ^r+r", they almost entirely cancel in the cross section times branching ratios 
and the remaining part is so small that it has no practical impact whatever benchmark scenario is considered. 
This can be seen from the almost identical tables XI-XIV and Figs. 4 of Ref. [14] that describe four benchmark 
scenarios [2T] . 
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a cross section that is only ^ 10% higher; this can be attributed to the different central scale 
and renormalization scheme for Mb that have been usecH. 
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Figure 1: The normalization of the cross sections ct^^a (left) and o"^^^ (right) at Tevatron 
energies as a function of Ma when using the MSTW PDFs and unit Abb couplings. In the 
inserts, shown are the various sources of theoretical uncertainties when the rates are normalized 
to the central values. 

For the evaluation of the theoretical uncertainties that affect the Higgs production cross 
sections as well as the decay branching ratios, we will proceed as follows. 

The Higgs decays branching ratios (BR) have been discussed in Ref. [15] and are simply 
affected by the parametric uncertainties on the input 6-quark mass and as- As the QCD 
corrections to the dominant ^ ^ bb decays are large, they are resummed by switching from 
the 6-quark pole mass M^ which appears at tree-level to the running quark mass in the MS 
scheme evaluated at the scale of the Higgs mass, m6(M$). The uncertainties come then from: 
i) the starting point for the 6-quark mass fni,(jni,) = 4.19lQ'Jg GeV [23] where the central 
value corresponds to a pole mass of Mh = 4.71 GeV; and ii) the error on the QCD coupling 
as(M|) =0.1171 ±0.0014 at NNLO (the value adopted in the cross sections) [H] which is used 
to run the 6-quark mass from mi, up to M<i,. Assuming that there is no uncertainty in the 
$ — )■ T~^T~ decay as only electroweak effects are present, and adding the errors in quadrature, 
one finds an uncertainty of ^ +4%, -9% on BR($ ^r+r") and ^ +1%, -0.5% on BR($->66) 
at the la level over the entire relevant Higgs mass range. Ma = 90-200 GeV. 

In the case of the production cross sections, the uncertainty from the missing higher orders 
in perturbation theory is usually estimated by varying the renormalization and factorization 
scales in the domains fio/n < fiR,fiF < '^/^o around the central scales /io, with the additional 



^We thank S. Heinemeyer for a discussion on these points. 
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restriction 1/k < fin/fip < imposed. While we choose k = 2 for the gg ^ A process, the 
value K = 3 is adopted for bb A. The reason is that it is well known that the cross sections 
in bb ^ A and in the twin process qq, gg — )■ bbA in a four-flavor scheme differ significantly [21] 
and by allowing for a wider domain for scale variation and, hence, a larger scale uncertainty, 
the two results become more consistent with each other. Furthermore, in the gg A process, 
there is an additional uncertainty that we will consider: the one due to the choice of the scheme 
for the renormalization of the 6-quark mass. The latter is estimated by taking the difference 
between the results obtained in the on-shell mass and MS schemes and allowing for both signs. 
The inclusion of this additional effect is similar to increasing the domain of scale variation from 

K = 2 to K = 3. 

For the combined uncertainties from the PDFs and the QCD coupling Og, we will use the 
scheme made available by the MSTW collaboration [18]. The PDF+A^^Pa^ uncertainty, with 
a,(M|) = 0.120 ± 0.002 at NLO for gg ^ A and a,(M|) = 0.1171 ±0.0014 at NNLO for 
bb A, is evaluated at the 90% CL. To this, we add in quadrature the impact of a theoretical 
error on a„ estimated by MSTW to be A*^a, ^ 0.003 at NLO and A*^ ^ 0.002 at NNLO. 
Finally, in the case of the bb ^ A process, there is an effect induced by the uncertainty in the 
value of Mb in the 6-quark density. This effect is again estimated within the MSTW scheme 
by allowing for an uncertainty on the pole 6-mass of ±0.25 GeV from the MSTW central 
value Mh = 4.75 ± 0.25 GeV. This resulting uncertainty is also added in quadrature to the 
PDF+A^^Ptts one. Note that we have evaluated both cross sections with four other PDF sets 
and found that the maximal values are obtained with MSTW in gg A, while some other 
schemes give ~ 20% lower rates. 

Finally, there is the uncertainty on the 6-quark mass that affects the amplitudes of these 
processes, in much the same way as what has been discussed for the Higgs decay branching 
ratios. It is estimated by evaluating the maximal values of the cross sections when one includes 
the error on the input MS 6-quark mass at the scale mf,, mf,(mb) = 4.191q;o6 GeV, and in the 
case of the bb ^ A process where the Yukawa coupling is defined at the high scale, cx fnh{fiR), 
the error on the coupling q;s(M|) = 0.1171±0.0014 at NNLO (in this case, it is almost the same 
uncertainty as in the A ^ bb decay). 

The results for these uncertainties on the production cross sections at the Tevatron are 
shown in Fig. [T]for the gg^A and bb^A processes as a function of M^. In the gg case and 
almost independently of Ma, the scale variation in a domain with k = 2 leads to an uncertainty 
(9(±20%), while the uncertainty from the scheme dependence in the renormalization of Mb is 
about ±6%; they add up to ~ 25% that is only slightly lower than the scale uncertainty in 
the bb process, ~ 30% for low Ma, in which the domain of variation is extended to k = 3. 
In the gg A (bb A) channel, the POF+A'^^P+^'^a, (with AM;, in addition for bb A) 
uncertainties are at the level of ±10% (±20%) for Ma ~ 100 GeV and larger (±30% in bb-^A.) 
at Ma ~ 200 GeV where the more uncertain high Bjorken-x values for the gluon and bottom 
quark densities are probed. The parametric error on leads to a ~ ±13%, —4% uncertainty 
in the gg —> A process and slightly less in the case of bb A. 

We turn now to the issue of combining these uncertainties. Clearly, the scale and scheme 
uncertainties, which are purely theoretical and both emerge from the truncation of the pertur- 
bative series, should be added linearly. The PDF+Aas+Arrib uncertainty, that we also would 
like to view as a reflection of the theoretical ambiguities due to the parametrization of the 
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PDFs and which have no statistical ground, will be evaluated on the minimal and maximal 
values of the cross sections with respect to scale and scheme variation (see Ref. [15] for the 
argumentation). This procedure gives results that are similar to those obtained with a linear 
addition of the scale+scheme and PDF+Aa^+Amb uncertainties as advocated in, for instance, 
Ref. [25] • Finally, at this stage, we add linearly the parametric uncertainty on mb which, in 
the case of interest, will drop anyway in the final result (see below). 

The combined uncertainties on the cross sections, when using this procedure, are also 
shown in Fig. [1] for the two production channels. As can be seen, they are very large: at low 
Higgs masses. Ma ~ 100 GeV, one has +55%, -35% for a{gg A) and +60%, -40% 
for a{bb — )■ A) which become at masses Ma ~ 200 GeV, respectively, ~ +60%, —40% and 
^ +50%, -40%. 

To consider the final state topology that has been searched for by the DO and CDF col- 
laborations [H], i.e. pp —)■ Higgs —>■ r"''r~, one has first to add the cross sections for the two 
channels gg ^ A and bb — )■ A, and then to multiply the resulting production cross section by 
the Higgs branching ratio BR(yl — t- r"^r^) ~ 10%. The resulting a{pp^ A) x BR(/1 -> r+r") 
at the Tevatron is shown in Fig. [2] as a function of Ma- We stress again that to obtain the 
true rate for ^ = A + H{h), one has to multiply the given values by a factor 2tan^ /3. 

In Fig. 121 shown also are the associated overall theoretical uncertainties. The uncertainty 
from the cross section alone is dominated by that of gg ^ A at low Higgs masses and bb—^ A 
at high masses as the corresponding cross sections are largest. In the product a{pp A) x 
BR(y4 T^T~), the parametric uncertainty that its common to the production and decay 
rates almost cancels out as shown by the solid curves in Fig. |2]and only a few percent are left. 
This leads to a smaller uncertainty in a{pp^A) x BR(y4 — )• r+r~) than in a{pp^A) alone. 
The final theoretical uncertainty for pp^ A^t^t~ at the Tevatron is of order +50%, —40%. 

To illustrate the impact of these theoretical uncertainties on the MSSM [M^,tan/3] param- 
eter space that is probed when searching experimentally for the pp^^ ^ t^t^ channel, we 
show in Fig. |3]the contour of the cross section times branching ratio in this plane, together with 
the contours when the uncertainties are included. We apply the model independent 95%CL 
expected and observed limits from the CDF/DO analysis (Table X of Ref. [H]). However, 
rather than applying the limits on the central ax BR rate, we apply them on the minimal one 
when the theory uncertainty is included. Indeed, since the latter has a flat prior, the minimal 
CTX BR value is as respectable and likely as the central value. One observes then that only 
values tan /3 > 50 are excluded in the mass ranges, M$ 95-125 GeV and M$ > 165 GeV. In 
the intermediate range M$ ^ 125-165 GeV, the exclusion limit is tan /3> 40-45, to be con- 
trasted with the values tan/3>30 excluded in the CDF/DO analysis. Hence, the inclusion of 
the theory uncertainties has a drastic impact on the allowed [MA,tan/3] parameter space. 

Finally, let us note that there is a subleading channel which has also been considered at 
the Tevatron, fe(yf — )■ $6 — )■ 3b [26]. The evaluation of the theory uncertainties in bg ^ ^b [9] 
follows that of the parent process 66 — )■ $ (for which it is part of the NLO contributions) and 
similar results, i.e. a total uncertainty of ~ ±40%, are expected. However, in this case, it is 
the ^ ^ bb decay which is considered experimentally and, since BR($— 7-66) has a small error, 
the uncertainties in a and a x BR are almost the same. These uncertainties will thus also 
impact the excluded [MA,tan/3] parameter spac^. 

^Notc that the 2a excess observed by CDF in this channel cannot be a Higgs signal as it would correspond 
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Figure 2: a{pp ^ A) x BR(A r^r^) as a function of Ma at the Tevatron, together with 
the associated overall theoretical uncertainty; the uncertainty when excluding that on the 
branching ratio is also displayed. In the inserts, shown are the relative deviations from the 
central values. 
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Figure 3: Contours for the expected a{pp^^ ^t^t ) rate at the Tevatron in the [M^,tan/3] 
plane with the associated theory uncertainties, confronted to the 95% CL exclusion limit. 



to a much larger excess in pp — > r+r^ which has not been observed. In addition, this 3b channel probes 
huge tan/?> 100 values, whereas values tan/3>50 are^heoretically not favored as they lead to a $55 Yukawa 
coupling that is non perturbative. For such values, the total Higgs widths are to be included. 



In conclusion, we have updated the cross sections for the production of the MSSM CP-odd 
hke Higgs bosons $ at the Tevatron in the processes gg ^ ^ and bb ^ ^ and found smaller 
rates in the high Higgs mass range compared to those assumed by the CDF and DO experiments. 
We have then evaluated the associated theoretical uncertainties, including also the ones in the 
$ r^r^ decay branching fractions, and find that they are very large. These uncertainties, 
together with the correct normalization, affect significantly the exclusion limits set on the 
MSSM parameter space from the negative Higgs searches in the channel — $ — )■ r~^r~ at 
the Tevatron. Additional material is given in Tables 1 and 2 which summarize our results. 
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1 n <^'7 +3.8% 

iU.O/ _g.^^^ 


n 1 +51.8% 
^•^^ -38.5% 


1 fin 

iDU 


1 /in 

i.4U 


+21.2% +5.4% +12.9% +13.4% +57.1% 
-19.2% -5.4% -13.2% -4.1% -38.9% 


1 n fi'? +3.8% 
10.63 _g75,^ 


n 1 c: +52.0% 
^■^•^ -38.9% 


iDO 


11/1 
1.14 


+21.2% +5.4% +13.1% +13.4% +57.4% 
-19.2% -5.4% -13.4% -4.1% -39.1% 


1 n fi8 +3.9% 

10.68 _g75,^ 


n 1 9 +52.3% 
^•^^ -39.1% 


1 7n 




+21.2% +5.4% +13.3% +13.4% +57.6% 
-19.3% -5.4% -13.6% -4.1% -39.3% 


1 n 7A +3.8% 


n. , r. +52.4% 

^■^^ -39.4% 


175 


0.76 


+21.2% +5.4% +13.5% +13.4% +57.8% 
-19.3% -5.4% -13.9% -4.2% -39.6% 


10.79 tlS^ 




180 


0.63 


+21.1% +5.4% +13.7% +13.4% +58.1% 
-19.4% -5.4% -14.1% -4.1% -39.8% 


10.84 +3.95 


n 07 +53.1% 
-39.7% 


185 


0.52 


+21.1% +5.4% +14.0% +13.4% +58.3% 
-19.4% -5.4% -14.3% -4.2% -40.0% 


10.90 tl'S. 


06 +53.1% 


190 


0.43 


+21.1% +5.4% +14.3% +13.5% +58.5% 
-19.5% -5.4% -14.5% -4.2% -40.3% 


10.95 tl^^ 


n r\K +53.3% 

U-UO _4o.4% 


195 


0.36 


+21.1% +5.4% +14.6% +13.4% +58.8% 
-19.5% -5.4% -14.7% -4.2% -40.5% 


1 1 00 +3-8% 
11. UU _ggy^ 


'^•U* -40.6% 


200 


0.30 


+21.1% +5.4% +14.9% +13.4% +59.0% 
-19.6% -5.4% -15.0% -4.2% -40.7% 


1 1 04 +3-9% 
11. U4 _g.7^^ 


n nq +54.2% 
U-UO _40.7% 



Table 1: The production cross sections in the gg ^ A process at the Tevatron (in fb) for 
given A masses (in GeV) at a scale fip = fiji = \Ma with MSTW PDFs. Shown also are the 
corresponding uncertainties from the various sources discussed as well as the total uncertainty. 
In the other columns, displayed are the branching ratio BR(y4 — t- r+r^) [in %] and the product 
(j{gg yl)xBR(A — r'^r^) together with their respective total uncertainties. 
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Ma 


NNLO 
"bb->A 


scale PDFs param total 


BR(tt) 


o-xBR 


90 


26.31 


+35.1% +18.2% +10.5% +70.4% 
-29.9% -17.2% -3.3% -44.9% 


Q 60 +3-8% 


9 c-o +63.4% 
^■^•^ -44.0% 


95 


21.04 


+30.6% +19.0% +9.9% +64.1% 
-29.0% -16.6% -3.7% -44.9% 


Q 70 +3-7% 
^•'^ -8.8% 


^•^^ -43.9% 


100 


16.96 


+26.6% +18.7% +10.2% +61.5% 
-27.0% -17.5% -3.4% -42.9% 


q 70 +3.8% 
^•'^ -8.8% 


1 r.. +54.5% 
-42.0% 


105 


13.78 


+22.6% +19.2% +10.7% +57.8% 
-27.2% -17.6% -3.4% -43.4% 


Q +3.8% 


1 oc +50.8% 
-42.4% 


110 


11.22 


+20.3% +19.4% +10.1% +54.7% 
-25.5% -17.9% -3.6% -42.6% 


Q 96 +3.8% 


1 19 +47.9% 
^•-^^ -41.7% 


115 


9.18 


+ 17.6% +20.1% +10.6% +52.7% 
-24.8% -17.9% -3.6% -41.6% 


1 n 04 +3.8% 
iU.U4 _gg^^ 


Q 09 +45.7% 
^•^^ -40.6% 


120 


7.57 


+ 15.1% +21.0% +10.9% +50.2% 
-24.1% -17.8% -3.0% -40.3% 


10 19 +3.8% 


n 77 +43.1% 

^- ' ' -40.0% 


125 


6.29 


+ 13.6% +21.5% +10.2% +47.9% 
-23.5% -18.1% -3.2% -40.4% 


10 1Q +3-8% 


0-64 ^ 


130 


5.24 


+ 12.8% +21.7% +10.5% +48.0% 
-22.9% -18.6% -3.2% -40.4% 


1 96 +3-8% 




135 


4.36 


+ 12.4% +21.8% +10.5% +48.5% 
-21.9% -19.0% -3.3% -39.9% 


1 n +3.8% 


0-45 „38:9%o 


140 


3.66 


+ 11.6% +22.5% +10.2% +46.9% 
-21.5% -19.1% -3.6% -40.2% 


1 3Q +3-8% 


n 00 +39.8% 
-39.3% 


145 


3.08 


+11.0% +23.0% +10.4% +46.7% 
-21.2% -19.5% -3.3% -40.3% 


10 46 +3-7% 


n 09 +39.5% 
^•"^^ -39.4% 


150 


2.60 


+ 10.8% +24.1% +10.4% +47.3% 
-20.3% -19.3% -3.4% -39.7% 


1 ^9 +3.7% 


n 97 +39.9% 
-38.9% 


155 


2.20 


+ 10.8% +24.3% +10.3% +47.8% 
-20.0% -20.1% -3.5% -39.3% 


1 ^7 +3.8% 
iU.O^ -8.7% 


n 90 +40.8% 
'^•^■^ -38.2% 


160 


1.88 


+ 10.5% +24.9% +10.4% +47.8% 
-20.0% -20.4% -3.4% -40.0% 


10.63 tl^^ 


0-20 


165 


1.60 


+ 10.6% +25.1% +10.4% +49.2% 
-18.7% -20.6% -3.3% -38.5% 


10.68 +3.9| 


n 1 7 +41.8% 

-2,7 AVo 


170 


1.37 


+ 10.6% +25.3% +10.5% +49.8% 
-18.4% -21.1% -3.1% -38.7% 


1 74 +3-8% 


p. 1 c +42.4% 
^•^^ -37.8% 


175 


1.17 


+ 10.5% +26.0% +10.4% +50.5% 
-18.2% -21.5% -3.4% -38.8% 


10.79 +3.8J: 


n 1 q +43.0% 
^■^■^ -37.6% 


180 


1.01 


+ 10.3% +26.7% +10.4% +50.7% 
-18.1% -21.6% -3.4% -38.6% 


10.84 


n 1 1 +43.2% 
^•^^ -37.4% 


185 


0.87 


+ 10.4% +27.8% +10.3% +51.4% 
-17.6% -21.5% -3.2% -38.6% 


10.90 tlS^ 


OQ +43.8% 
^•^^ -37.6% 


190 


0.75 


+ 10.3% +27.2% +10.4% +52.1% 
-17.3% -22.7% -3.3% -39.2% 


10.95 tlS^ 


n ns +44.4% 

-38.2% 


195 


0.65 


+ 10.5% +28.5% +10.1% +52.3% 
-16.2% -22.6% -3.4% -38.9% 


1 1 00 +3-8% 




200 


0.56 


+ 10.4% +28.5% +10.6% +53.2% 
-16.3% -23.6% -3.2% -38.9% 


1 1 04 +3-9% 


0.06 



Table 2: The production cross sections in the bb ^ A process at the Tevatron (in fb) for given 
A masses (in GeV) at a scale Hf = I^r = with MSTW PDFs. Shown also are the 

corresponding uncertainties from the various sources discussed as well as the total uncertainty. 
In the other columns, displayed are the branching ratio BR(y4 t^t~) [in %] and the product 
a{bb — > ^)xBR(A — )■ t^t~) together with their respective total uncertainties. 
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